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ABSTRACT 

To  solves  Fy  = 0 numerically  we  use  two  different  methods,  the  first 

^ t ^ ' . f 

of  which  is  sketched  already  Secondly  introduce' a neighbouring 

problem  (N.P.)  Fu  - d l|d  ||  small,  with  known  solution.  ^We  solve  the 

i-  * . i 

original  problem  and  the  N.P.  with  theTsame  discretization  method"'.'  The 
known  error  of  the  N.P.  is  used  as  an  estimation  for  the  unknown  error  of  the 
original  problem.  These  procedures  are  used  iteratively  and  their  relations 

^ ^ • «-v 

are  discussed.  In  subsequent  papers  we-wiH  apply  our  general  theory ^to 
some  special  cases  and  will  discuss  relations  to  collocation  metliods  and 
to  PEREYRA's  deferred  correction  methods, [ 8,9]  . 
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Defect  Corr'ections  via  Neighbouring  Problems 
1. 6encx~-al  Theory. 


Klaus  Bohmer 


Summary : To  solve  a functional  equation  Fy  = 0 numerically  we 
use  two  different  methods:  The  first  one,scetched  already  in 
[3],  combines  Newton’s  method  with  a discretization  of  the 
linear  problem.  Secondly,  we  intr>oduce  a neighbouring  problem 
(N.P.)  Fu  = d,  II  d 6 small,  with  known  solution.  We  solve  the 
original  problem  and  the  N.P.  with  the  "same  discretization 
method".  The  known  error  of  the  N.P.  is  used  as  an  estimation 
for  the  unknown  error  of  the  original  problem.  Both  methods  are 
used  iteratively  and  their  relations  are  discussed.  The  idea  of 
using  N.P.  goes  back  to  ZADUNAISKY  [15,16,17],  again  discussed 
by  STETTL’R  [14].  They  treated  initial  value  problems  of  ordinary 
differential  equations  where  STETTER  used  our  first  method,  too. 
The  method  of  N.P.  was  applied  by  FRANK,  HERTLING,  UEBERHUBER 
[5,6,7]  to  initial  and  boundary  value  problems  of  ordinary 
differential  equations ., In  subsequent  papers  we  will  apply  our 
general  theory  to  some  special  cases  and  will  discuss  relations 
to  collocation  methods  and  to  PEREYRA’s  deferred  corrections 
[10,11].  We  find  numerical  methods  which  seem  to  be  as 
convenient  and  appropriate  than  those  derived  directly  (see 
[5,6,7,14,15,16,17] ). 


1.  Introduction  " 

Witli  Stetters  [13]  notations , slightly  modified, a discretization 
method  applicable  to  the  given  problem 

(1.1)  F(y)  = 0,  F : D •>  E°,  D E,E°  Banach  spaces. 


is  an  infinite  sequence  ( ° »‘^’h^h  6 I]  ’ ^^^^0,1)^]  C 3R  + , 
inf  IH  =0,  such  that 


1 


V 4 • « / 


; E -*■  E^,  A^°:  E°-f  linear  bounded, 

r.  Banach  spaces,  dim  R = dim  E ° < ~ 

li  ’ h n . h 

r is  continuously  imbeded  into  a Banach  space  E and 

limljAj^yll  = ||y||g  for  every  fixed  y € E S e 
h-*-o 


liml|Av,z|l  = ||z||  for  every  fixed  z e E . 

It+o 

*^h  • S ^^h  ^h°^’  C,  (E  ->  E°)  and  F 6 for  h e H 


More  !1‘ 11  means  anyone  of  the  norms  in  E,  E°,  E^^ , Ej^°. 

A'.'plying  jyi  to  problem  (1.1)  we  find  the  discretization  of  (1.1) 


(1.3) 


«h  ^ “h"  ’ “hS  Eh  for.  hen. 


»h(nh>  = 0. 


Wo  assume  that  (1.1)  and  (1.3)  have  unique  solutions  y € d 
and  6 Dj^  (see  [13]  pp  12  ff.).  Further  we  use  the  equivalent 


notations 


Xi  = X2  + 0(h’^),  r e ]R  ^ iff  ||x^  - X2II  = O(h^)  for  h-*-0, 


If  there  is  a sequence  of  Aj^  : E ->■  E , h 6 H such  that 
^h^^h'^^  ” ^h  ^F(u)  + A^(u)}  for  every  u 6 E 

Ajj  is  called  a "local  error  mapping".  In  many  cases  the  local 
error  mapping  admits  an  asymptotic  expansion  up  to  the 
order  Vq,  that  is 


( 1 . S ) 


A}^°Aj^(u)  = f h 'f^(u)  + 0(h  ^^^)},0<v^ 


<V2<. . 


for>  u 6 Dq  ^ D,  f^  :D^->E  ,f^  independent  of  h. 


2 


If  y e.  Dq  is  the  solution  of  (1.1)  and  if  (1.3)  is  consistent 

to  (1.1)  of  order’  v then  the  "local  discretization  ei-'ror" 
o ^ 

:=  satisfies 

(1.6)  = 4>^(Aj^y)=A^°{  f h^f  iy)  + 0(h'^‘^^^)>  . 

i =p 

Very  important  for  numerical  applications  is  the  question, 
whether  (1.6)  carries  over  to  tire  "global  discretization 


(!•?)  Yj^  :=  hj-j  - ^}^yj  y solutions  of  (1.1),  (1.3). 

Gragg  [8]  has  studied  this  question  first  for  initial  value 
problems  of  ordinary  differential  equations.  Stetter  [12,13] 
general:! zed  Gragg's  result  to  functional  equations.  In  these 
papers  we  always  have  = r or  = 2i.  Several  difficult 
special  cases  were  treated  directly  (see  for  example  Benson 
[1]  ). 

2.  Asymptotic  expansion  of  the  global  discretization  error 

Here  we  generalize  and  modify  Stetters  [12,13]  results  a little 
bit  to  prepare  it  for  our  later  applications.  We  require  that 
the  in  (1.5)  are  such  that 


(2.1) 


IN..  - ’•  • • ’'’q»''q+l^  ’ • ^^q^^q  + 1 


q q+1 


V i ,v  j e I'l  5c  implies  ^ orv^  ^ '^q+l 


for  this  case  we  give  (see  Stetter  [13],  p.  25). 

Definition  2.1:  We  call  the  asymptotic  expansion  (1.5)  of  the 
local  evror  mapping  (,v^^v^)-smooth  at  u,  if  the  derivatives 

(a)  V -V  V -V 

f^  (u ) ,q  = 1 , 2 , . . . , [—2 — ^]=max{n  £ Z|n  < — ^ exist  and  if 

for  u,e.  £ D =D  C.  D,||  ^ h ^ e |]  small  enough^  k = p(l)q,  the 

^ q ^=P  ^ 

following  relation  is  valid 


f 


m*  T 


( ? . 2 ) 


V -V, 


^=1  a=i 


(u)(  i h 6j^)  } 
k=p 


+ ? h'^^e  ) + 0(h  . 


= t '•  ^V: 

^ k = p ^ 


In  many  cases  the  conditions  u,Cj^  6 D = D are  too  stringent 

4 V 

for  (2.2).  So  we  give 


A 

i 

1 

> 


1 


i 

\ 

1 

■j 


Definition  2.2:  Let  D„  ,,  c k = p(l)q,  be  such  that 

' — ■ 4>^ 

D^D  ,£D,D  £E,  and  that 

q q,k  ’ q q’ 


A,  (u  + ) h ''  e,  ) and 

^ k=p  ^ 


(2.3) 


fj^^u)  i I ^ y x=l(l)q,  o = Hl)[^^— 


q"  V 


k = p 


with  q”  = q"(q}p,i,o)  suoh  that  v^.;  = (o-l  )Vp 


are  defined  for  u 6 D , e,  6 . Further  let 

q K q , 


(2.4)  S 


^j^'*(u  + f h ^ e.)  :=  f h ^(f,  (u)  + 
k = p ^ 


\=1 


V -V 

[-9--1] 

1 f(0).  ..  , 

+ I ^ 

0 = 1 


k=p 


. k nO^ 

h o,^)  ) 


eatiefy 


(^.S)  Aj^(u  + ? h c^)  = 


9 V. 


k = p 


j^/  - Aj^''(u  + I h e.  ) + 0(h  9 
k = p 


- ^ - 
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are  called  ( aoymptoticallij ) admissible  sets 


Then  the  D , 

q»k 

(for  Wl  and  ^ }. 


The  summation  index  q”  in  (2.4)  is  chosen  such  that  with 

^ G + 1 V 

k > q-'  only  contributes  to  0(h  ^ ).  Since  fui^ther  < 


(a) 


P = Q” 


iff  V < V -ov  , only  those  f (u)  contribute  to  A,"  with 
l=qp’-'  1 n 


V,  < v-av^. 
i - q p 


Now  the  following  Proposition  is  simply  proved: 

Proposition  2.3:  Let  the  operators  given  in  (2.3)  be  defined 

and  continuous  for  u 6 D , e,  6 D , with  respect  to  the  norms 

q K q 

in  E , E , . Then  (2.5)  is  satisfied. 

q q>k 


Since  we  are  concerned  only  with  asymptotic  expansions, 
h ^"'"^-terms  are  fairly  uninter'esting . So  in  addition  to 
ordinary  Frechet-derivatives  we  introduce  the  following  type : 

Definition  2.4;  is  called  v-times  ,V[^)-diffei'‘entiabte 
at  ri,  if  there  are  bounded  a-linear  operators  (n)-' 


(EjP  E,°,  such  that  for  V >rv  , e 6 D , \ -p(l  )q,  cnc? 
nn  HP\q>i 

q 

II  ^ h e 11  small  enough 
\=p  ^ 


(2.6)  < 


4>h  (n  + I h 'e^)-  Oj^(n) 

I =p 


0=1  ^ ^ 1=P  ' 


m3 


0(h 


in{ 


The  operator  called  the  o-th  (v^^\'^)~derivative  of 

at  n • 

Now  the  following  Proposition  is  a simple  consequence  of  the 
last  two  definitions: 


- 5 - 


Proposition  2.5;  Let  smooth  in  B^Cy)  ; ={u  € Dj  ||u-y l|<p} 

and  lett  with  v > F have  r Lipsahitz-continuous  derivatives . 

Then  {F+Aj^}  is  r-times  (.v^,v^)~differentiable  in 

and  with  j such  that  v ^ <min{ (r+1 )Vp -aVp  we  find 

(2.7)  (Aj^u)  = A^°{F^°\u)  + J h ^ f|”\u)},  a = l(l)r. 

For  o < r < with  at  most  one  =~sign,  we  always  have 

min{ (r+l)Vp,v^^^}-aVp>Vp,  so  j > p > 1 and  the  I in  (2.7) 
is  well  defined.  .For  a = r = that  j < 1. 

Then  we  define  I :=  0 in  (2.7). 

The  proof  of  Proposition  2.5  is  obvious  with  Def inition;  2 . 1 
and  2 . 4 and  (2.1). 


To  give  the  following  Theorem  we  need  some  formal  notations. 
Let  F be  [ v^/v^] -times  differentiable,  let  y be  the  solution 
of  (1.1)  and  let  A^  be  (v^,Vp)-smooth.  Then,  with  F(y)=0, 
f^(y)=0,  i=l(l)p-l,  y,  e^e  D 


q V .i  9 

(F  + Aj^)(y  +1)1  e^)  + 0(h  ^■^^)  = 1 h 


(2.8) 


[v,./v  ] 


V <v  -av. 


+ I I )i  e ) 

rt-i  •*-'1  ^ ^ 


Now  we  introduce  operators  g j (.,...,. ) by  collecting  equal 
powers  of  h in 


V =v  ~0\L 
^ q P V 


k = p 


(2.9) 


q V: 


= I b ^g^(6p,. . . ,e^^jj),  gj  = 0,  j=p(l)p^-l. 


with  V “ M 'i  " '’-i  ” » )^bercfore  S,(j)  < j; 

P J P 


6 


"it  I I - - - - ^.■...  , ...  ■■  l.-VJI. 


in  the  first  line  of  (2.9)  ej^  gives  contributions  to  powers 
> Vp  + Vj^.  So  the  biggest  k appearing  in  g^  (ep , . . . ,6j^)  is 

k=JUi)  with  V.  = V + V...V.  For  the  special  case  v =\,  i=l(l)q, 
..  3 p JU3)  ^ \ 

we  have  ^-(j)  = j - p.  With  the  g^  , defined  in  (2.9),  we  finally 
have 


(2.10)  i 


(F  + A^)(y 


o(h  = 


? .“kF 


I =p 


:=p 


' (y)6j+f j (y) 


V,  +V,  =V  . 
I k j 


’•  • • ’®)l(j)^^  ’ ^)Uj) 


Since  we  want  to  have 


n = Aj,(y 


F f 


I =p 


V V . 

h ^ e^  ) + 0(h  ^■^^) 


and  since 


<&j^(Aj^u)  = A^°{F(u)  + A^(u)} 

the  proof  of  the  following  Theorem  mainly  consists  in  equating  to 

zero  all  the  coefficients  of  h , i=p(l)q,  and  can  be  done  exactly 

in  the  same  way  like  in  Stetter  [13J.  (Due  toll,  v ([  v /v  l+l)>v 

" p q p = q+1 

Theorem  2.6:  Let  the  original  problem  ^ with  the  exact  solution  y 
and  the  discretization  method  111  , applicable  to  ^ ^ satisfy 

(a)  ^ is  stable  for  ; (see  [13]); 

(b)  M is  consistent  with  of  order  v ; 

(/  p 

(c)  an  asymptotic  expansion  (1.5)  of  the  local  error  mapping 

A.  exists  and  y 6 D ; 
h Q 


(2.11)'^ 


(d)  Aj^  is  -smooth  at  y; 

(e)  F has  [v^/Vp]  Lipschitz-continuous  hounded  derivatives 

in  B^(y); 

P 

_ 

(f)  (F'(y))  exists  and  is  bounded. 


7 - 


Define  the  , j = p(l)q  by  (see  (2.9),  (2.10)) 


(2.12)  F'(y)e.  = -{f-(y)  + I f ' (y)  e.  + g.(e  , 

^ \^l,k|p  ^ ^ J P 

V^PVk=Vj 


and  let  e-  € . 

3 q>3 


Then  the  global  discretization  error  admits  a unique 

asymptotic  expansion 


q V +-1 

(2.13)  n,  - A^y  = A^  I h^e  +0(h‘^^) 

\=p 


For  the  rest  of  the  paper  we  postulate  (2.13) 


For  our  later  applications  we  need  the  structure  and  the  arguments 
of  the  gj  in  (2.9)  a little  bit  more  completely: 

Proposition  2.7:  The  gj  in  (2.9)  are  operators  acting  multi- 

linearly  on  the  e,  , k = p ( 1 ) it  ( j ) , v-  = + Vp/--.,  and  linearly 

on  the  o = 2(l)[v-/v  ],  and  the  f^°'^(y),  o = 2(l)[— ^ b]  . 

(,\  ^ ^ v.-v  ^ '^p 

f''  (y)  contributes  to  g-  iff  2 < [— ^ i-]  or  for  i = 2(l)\  (j) 

I 3 - Vp 

with  ^ v.-2Vp  < Sr 


([v./v  ]) 

gj=g4(e  ,.  . . ,e^(  .);F"(y)  ,.  . . ,F  ^ ^(y);f|[(y), 


(2.14) 


v.-Vi 


([  -\i  ) 


...,fj  P (y) ,f" (y) > • • • >f ^ 


.,f;'(j)(y)) 


Here  F^^^  (y)  I’esp.  f,^°\y)  has  to  be  applied  to  those  q with 

I K. 

(2.15)  V.  < v.-(a-l)v,^  resp,  v,  < v-  - v - (o-l)v  . 

N-J  P K-]l  p 
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Pi-'oof : We  only  give  the  arguments  for  the  By  (2.9)  we 

have  contributions  to  h ^ g.  by  \y)(  f h ^ 

k = p 

^ v_.-v, 

, o r J ^ ■ 


V +av  < V.,  that  is  for  a = 2(l)[-t^; . From  2 < [ 

I P = 1 ’ 


we  get  \ ( j ) . 


3.  Defect  corrections  via  Newton's  method 

The  results  given  here  are  slight  modifications  and  extensions  of 
those  given  in  [3].  We  use  an  operator  (see  for  exampile  [9]) 


(3.1)  Tg  : 


, 6 e ]R  Q :=  {teiR  1 t > 0} 


p,  y,  : = T~n, 

'h  -^h  6 h 


whicli  essentially  reproduces  the  asymptotic  expansion  of  h,  in  (2.13) 


(3.2) 


If  n,  satisfies  (2.13)  then 
h 


?v  V- 

h ^ (y)  + 0(h  *^)  , 


'^6'^h  = ^h  = y ^ ^ e^(y)  + 

^ =P„ 

where  v—  ^ v„~S  < < v_  . . 

q - q q = q + 1 


Usually  T.  are  interpolation  operators  and  often  1|*||^  includes 
c ^ 

derivatives,  so  one  has  q < q.  We  will  have  to  compute  the  defect 
F(y  ) = F(T  n,  ) and  we  need  an  asymptotic  expansion 


(3.3)  F(y,  ) = 


V V" 

h ^g  (y)  + O(h^) 


- 9 - 


with  the  q from  (3.2)  (see  (3.7)).  Sometimes,  for  example  in 
spline  interpolation,  is  not  smooth  enough  to  admit  the 

full  expansions  in  (3.2)  resp.  (3.3).  In  those  cases  we  have  to 
cut  down  the  asymptotic  expansion  (2.13)  to 


’^h  ~ ^h^  = ^h  ^ ^ ^ 

1 =p 


) with  Qq  < q 


such  that  this  expansion  admits  (3.2),  (3.3). 

As  a consequence  of  the  modified  Newton  method  we  find  [2] 

Theorem  3.1:  Let  6 D he  an  approximation  for  the  solution  y of 

(1.1)  such  that  F is  at  least  K :=  -times  differentiable  loith 

uniformly  Lipsahitz-continuous  derivatives  in  B (y„)  and  let 
o . P ° 

F“(y^)  € oe  (E,E^)  be  invertible  with 

l!(F"(y^))"^||  < p,  p-sup  ||F'(z)-F"(y  )||  < q < 1 

p o 


||F"(y^)"^F(y^)||  < p*(l-q)  . 


Then  the  modified  Newton  method 


(3.4) 


F>=(yo)(y,-y,l-l>  = - ’ ^ = 1,2,3,...  , 


with  V :=  v,^ , V + v = v , £ = 1,2,3 

IN  ^ TN  T-..  TN  * * * 


P’  Po- 


£-1  Po  P£ 


defines  a sequence  y.  converging  to  y,  lim||yn-y!!„  = 0.  If 

* £->-00  ^ “ lu 

further  F‘(y)  exists  and  is  bounded  and 

/n  « \ « t 

q+1 


(3.5) 


y_  = y + i h ^ e (y)  + 0(h  ),  e (y)  :=  e (y), 

, -r.  ‘ ‘ s'-*  1 


1=P 

F"(yo)  = F'(y)  + 


\=p 


/ V /N 

^ C (y)  + 0(h  "1^1),  V-  + 

I ’a 


C (y)  6 (E,E°),  C (y)  independent  of  h. 


then  y admits  an  asymptotic  expansi on  of  the  form 


- 10  - 


iifiiMMdttiaiil 


(3.6) 


q V +1 

Yj,  = y + i Ps,  ^ 

'=P£ 


resp.  y^  = y + 0(h  ^ ) for  > q 


If  P'^Cy^)  = P'(yQ)j  the  second  part  of  (3.5)  follows  from  the  Taylor- 


f ormula. 


If  y^  and  satisfy  (3.5)  then  we  have 


F”(y^)  . F'(y^)  + ? h""'  ^^(y)  + 0(h"''i^^), 


^^(y)  e £ (E,E°). 


Proof : The  convei->gence  of  t]ie  modified  Newton  metliod  under  the 
restrictions  on  y^  and  F"(y^)  follovjs  by  usual  conclusions  (see  [2]). 
(3.6)  is  proved  by  induction.  Let  yji__2  P®  °F  form  (3.6)  with  Pj^ 

and  e „ replaced  by  Pn_.,  and  e p_^.  Then  we  have  with  Taylor's 
formula  and  because  of  tlie  Lipscliitz-continuous  - we  suppress 

the  argument  y in  tlie.  asymptotic  expressions  - 


F(yj^_^)  - F(y  i I h ' e + 0(h  '^''^)) 


q V V . ; 

= F(y)  4 (F'"(y  ) 4 F'(y)  - F"(y^))  I h ^ e ._i+0(h  ^ •^), 
O O \ % )6  ” X \ 

j 

^ F^j\y)  ^ i I 

3 = 2 i=p„  . ^ 


Since  (k  4l)Vp  ^ (K'4l)Vp  ::  |y  Vp  4 Vp  > and  v-jith  F(y)  = 0 


P£-l 

and  (3.5)  we  find 


P£"^ 


^3  V, 


(3.7)  F(y^_^)  --  P=(y^)  I h ^ + I h ^ 


^=P£-1 


with  e and  e „ independent  of  h.  Nov\>  (3.4),  (3.5)  and  (3.7) 

ijX>*'X  TjX 

i mply 


5 


11 


1 

•; 

I 


I 


f 


i 


i 

1 


s 

'if 

I 

t 

1 

i 


I 


j 

' 


I 

3 

I 

1 

3 

I 

•J 

,t 

] 

i 

1 

i 

I 


y = y * I h'''  (e  - I-(y„)-l  S ) * 0O.Vl,  . 

>=p, 

Ai^air.  (3.5)  implies 

-1  -1  ^ '*1  '’n  + 1 

F“(y  ) = F'(y)  ^ (id^o  + I h ^ C"  + 0(h  ^ 

with  C”  independent  of  h,and  a combination  of  the  two  last  results 
pi'cvos  (3.6).  The  last  statement  again  follov.;s  with  the  Taylor  for- 
mula. □ 


Since,  again,  (3.4)  is  usually  not  solvable  exactly,  wc  have  to 
use  a discretization  method.  The  crucial  question  is,  if  this 
discretization  method  reproduces  the  asymptotic  expansion  of  the 
y^^ . To  guarantee  that  we  need 

Definition  3.2:  In  (1.1)  let  Fu=F^u+d  reap,  Fu=F^u+F2U  with 

F^jF^  €»<^’(E,E°)  and  d independent  of  u.  Let  Itt  , applied  to  , give 


(3.8)  i 


L^(FiU^F^u)=$^(F^)A^u.^^(F2)A^u  . 

Then  we  call  TUft  a linear  discretization  method  for  ^ and  4>j^(F^  )A^u+Aj^d 
= 0 reap.  (F^  )Aj^u+?>j^  (F2  )A^u  = 0 linear  discretizations  of  F^u+d  = 0 resp. 

F^u+F2U=0. 
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Let 


(3.9) 


^il-1  ■"  ’ «-  = l,2,...  and 


be  a linear  discretization  of  (3.4).  Sometimes  we  use  F” (y^) =f ' (y^) j 
then  often  <}>yj"  is  the  Frechet  derivative  of  from  (1.3),  but  that 
is  not  true  in  every  case.  We  want  to  prove  a result  for  the  ^ 
in  (3.9)  similar  to  that  for  the  y^  in  (3.6).  Since  we  have  to  use 
in  (3.9)  the  not  only  the  lower  index  p^,but  also  the  upper 

index  q has  to  be  changed  for  every  step  corresponding  to  , and  we 
define  (see  (3.2)  and  (3.4)) 


L. 


4 

i 

i 

,'j 

4 

4 


i 


(3.10) 


V.  < V 


- 6 < V " < V 


£-1 


q,-Hl  • 


Now  the  following  Theorem  is  valid 

Theorem  3.3:  In  addition  to  the  conditions  in  Theorem  3. 1 let 
$j^'KAj^y^) Aj^u  = Aj^°d.  be  a stable  linear  discretization  for  F"(y^)u=d. 
Further  let  in  (3.9)  the  and  therefore  the  corresponding  defect 

V ^ Vi 

' % — r\  * 


(3.11) 


) , £=1,2,... 


\=p 


£-1 


P£-l  V. 


F(yj^_l)  = F"(yo){  I ^ ^ I h 'e^^j,(y)  + 0(h 

^"P£-l  ’ ■ '"P£  ’ ' 


be  such  that  (with  e^  £-1  j jj, 
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J 


V<Vo>s.t-i  = V 


implies 


e n ^ = A 

1 ,11-1 


C>  * ^ V*’ 

\ ,n-l  V . s „ . 


(3.12) 


with  V + V 


= v„  , V + V 


' ' S.t-1*!  = <1 


> V"  and 
*^^-1 


^ o'V  o,  j. 

V^Vo^S  = S,£  S=^h^S,£  ^ I ^ 

1 =p  ’ ’ 


+ 0(h  ^ ’^"^)  . 


Then,  for  h small  enough,  t\^  ^ from  (3.9)  satisfies 


V 


(3.13)  ’3b  £ = ^ ^ ^ ^)/or  p < q 

5 ^ — -n.  *■5 


resp. 


(3.14)  = A^y  + 0(h  ) for  p^  > . 

Proof : In  a totally  analogous  way  to  (3.7)  we  find  from  (3.2), 
(3.10)  , (3.11)  and 


yii-1 


V 

y + I ^ n-i^y^  ^ ^ 

^=Pn-i 


that 


(3.15)  P(yp_i) 


P£-l  V V ^ '’qV 

F"(y  ){  I h „ .(y)>  + I h „(y)+0(h  ) 

'=P£-1  ’ '"P£ 


With  the  stability,  (3.9)  and  (3.12)  we  conclude 


or  finally  with  (3.11) 


‘A  V 


'h,£ 


i=p 


resp.  (3.14)  if  p^^  > q^^.  D 

To  find  (3.13)  or  (3.14)  one  has  to  know  T(yj^_^)  exactly.  In  many 
applications  that  is  not  the  case.  This  difficulty  is  overcome  in 
the  following 


Theorem  3.4:  Define  the  £ hy  h"  ^ “ '^h  o 


and 


(3.16)  v'Vo’<'ii;,o  -"h.ii-i’  = - > 1=1.2.... 


where 


'^Jl-l  V. 


(3.17)  ♦<n“_,.i)  = 1 h 'f,^,.,(y))l0(h’^-h.. 

Further  let  the  conditions  of  Theorem  3.  3 be  satisfied  with 

’^h,£-l ’^1  ,£-l  (3.11)5(3.12)  replaced  by  n=j!^  £-1  £-1 

a.s.o.  Then  (3.13)  resp  (3.14)  are  valid  for  instead  of  £. 

Proof : The  proposition  follows  by  comparing  (3.15)  and  (3.17).  □ 

The  question,  how  exactly  ^ in  ° should  be  known, 

may  be  treated  exactly  in  the  same  way  as  it  is  done  in  [10]  by 
introducing  the  concept  of  approximate  solutions. 

4.  Defect  corrections  via  neighbouring  problems 

If  the  discretization  (1.3)  of  the  original  problem  (1.1)  is  good 
and  if  li  is  small  enough  then  y^  :=  ^ (see  (3.1),  (3.2),  (3.9)) 

will  be  a good  approximation  for  the  solution  y of  (1.1)  and  the 
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1 


r 


defect  :=  be  small.  We  assume  that  the  neip.hbourinp; 

problems  (N.P. ) 

(4.1)  F(u)  = d , lld|l  < e N.P.  for  (1.1) 

and,  by  discretization  of  (4.1), 

(4.2)  ^h^^h^  ^ ^ ^o  ^ 


are  uniquely  solvable.  So  y^  is  the  unique  solution  of 


(4.3)  F(u)  = F(y^),  if  |lF(y^)||  < c^  . 


Though  we  know  the  exact  solution  y^  of  (4.3)  we  use  the  corres- 


j ponding  discretization  to  compute  ^ from 


(4.4) 


»h«h,o>  = ''h  '■‘yo>  • 


If  IlFCy^)!!  is  small  enough  tlie  known  error 


?h,o  - ‘hyo  = - ''h.c  * > 


see  (3.2),  (3.9),  should  be  a good  estimation  for  t)ie  unknown  error 
rij^  Q - Aj^y  of  the  original  problem  (1.1)  and  of  its  disci'-etization 

(1.3)  and 


(4.6) 


%,1  ==  ^h,o  ■ ^^h,o  - ^.^o^  = ’^h,o  - ^^,o  " ’^h,o' 


+ 0(h  ”) 


- 1C  - 


should  be  a better  approximation  than  ^ and  the  defect  F(y^) 
with  :=  T^rij^  should  be  smaller  than  F(y^).  ^ in  (4.6)  is 

usually  different  from  Pj^  in  (3.9).)  So 


F(u)  = F(y^)  . = Aj,"F(y^) 


are  closer  to  (1.1),  (1.3)  than  (4.3),  (4.4)  and  the  knov>?n  error 


«h.l  ' ‘hS'i  = ♦ 0(h  1) 


should  be  a better  estimation  for  the  error  ^ - A^^y  than  (4.5) 
and,  corresponding  to  (4.6),  we  define 

'’q" 

^h,2  •=  ^h,o  - ^?h,l  ■ ^h^l^  = '^h,o  ■ ^^h,l  “ ^h,1^^0^^^ 

This  process  may  be  used  to  generate  the  following  iteration  method 


^h^’^h,o^  = ° 


(4.7) 


^«,-l  '■  "^6  ’^h,t-l 


L '^h,«.  ^h,o  " ^^h,£-l  " ’^h,ii-l^ 


> ^=1,2,3, 


Since  we  hope  that  the  will  approach  the  exact  solution  y,  the 

defect  will  gro\<i  smaller  and  smaller,  so  t,-l  ^ ^ 

will  differ  less  and  less  from  ^)  = o.  Therefore  a very  good 
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X.MV  to  compute  from  (4.7)  is  to  use  a modified  Newton  method 

with  starting  value 

(ii.S)  . i = l,2,...j  Cq  = = O , 

and  hopefully  lim  C-  = . 

i-»-oo  n,x,  I 


As  a consequence  of  (2.6)  and  (3.3)  we  will  have  = 

I h + 0(h  ) with  appropriate  ^ 1*  Similarly 

i = p. 

as  in  3.,we  may  use  in  (4.8)  instead  of  ^’]!j(nj.|  any  (Hj^q)  with 


(n.  ) = I’llCnv.  ) + I (Pk  ) + 0(h"'l‘^^), 

h,o  h h,o  ^ i h,o 


(4.9) 


e cf  (Ej^,E^°)  , x=p(l)q,  v;;  + . 


q P q 


If  we  do  so  and  compute  only  the  first  step  in  (4.8)  we  find  the 
following  simplified  version  of  (4.7) 


^h^^h,o^  = ° 


(4.10) 


^£-1  "^6  '^h,£-l 

^h'^’^hjO^  ^^h,J!,-l  “ '^h,o^  " '^h 


1 = 1 


To  be  able  to  prove  the  following  Theorems  we  need 

Def  init.ion  ‘t . 1 : Let  m be  yV  )-smooth  for  . Furthei'  let^  for 


u e Bp(y ) , 


(4.11) 


4>j;(Aju)t  = d 


be  such  that  the  discretization  of 


(4.12) 


by  Ifft 


(4.13) 


F'  (u)v  = d 


()>j^(F' (u))t  = Aj^"^d 


satisfies 

q V +1 

(4.14)  (l)j^(F'(u))  = <^l'(A^u)  + I h ^ (J)  (Aj^u)  + 0(h  ^^^)  , 

i=p 

with  ({>^  e (E,^,Ej^°)  independent  of  h resp. 


(4.15) 


(F'(u))  = <^>h(A,u)  + 0(h  "^^^)  . 


'h""  " 'h'^h" 

Then  is  called  a differen liable  rasp,  a strongly  differentiable 

discretization  method  for  ^ . 

If, in  addition  to  (4. 9), we  have 

/s 

q V. 

I 

1 =^P 


Q V V 

(4.16)  F”(u)  = F'(u)  + y h ^ F (u)  + 0(h  ^'‘■^),  F € ojf  (E,E°) 


then  V'je  have  an  analogous  formula  to  (4.14)  with  (J)^(F’(u))  and  <?j!^(Aj^u) 
replaced  by  <{)j^(F”(u))  and  most  of  our  later  applications 
we  have  the  situation  given  in 

Definition  4.2:  Let  (Aj^u)f  n (4.9)  and  F”(u)  in  (4.16)  be  such  that, 
for  d and  u smooth  enough,  the  solutions  x and  v of 


(4.17) 


(Aj^u)x  = Aj^  d and  F”-(u)v  = d 


19  - 


satisfy 


?V  '^4.1 

h ^ f^(v)  + 0(h  ^ . 

i=p 

Then  F"(u))  is  called  ^ -admissible  ( 7^  a 

(.V  ,\)  )-smooth  discretization  method). 

It  is  possible,  to  give  sufficient  conditions,  depeiiding  on  d,u 
and  the  in  (4.9)  and  the  in  (4.16),  for  the  TfC  -admissibility 
of  (1>j^"(u),  F"(u)).  But  since  these  conditions  ai^e  very  complicated 
and  do  not  save  too  much  work  in  the  special  cases  to  be  treated 
later  on,  we  do  not  formulate  the  corresponding  theorem.  In  our 
later  applications  we  will  often  use 

4>j^KAhU)  = <t>jl,(A^u)  or  = *^h^' 

If  '7ft  is  strongly  differentiable  for  ^ , then  (^’^(Aj^u)  ,F' (u) ) 

is  "^-admissible. 

Theorem  4.3:  In  addition  to  the  conditions  of  Theorem  2.4  and 
3.1  let  ^ he  a linear  discretization  for'  ^ 

he  ~7dX -admissible  and  »Fj^°  (n^  o^^heH  stable  (see 

[13]  or  [10].),  that  is  for  this  case 

IIyII  i *^l^h'‘ ^'^h,o^'Yll  all  y e and  a fixed  c 6 IR^  . 

Further  let  (4.2)  be  the  dft -discretization  of  (4.1)  and  let  the 
N.P.s  (4.1),  (1.1)  resp.(4.2),  (1.3)  he  uniquely  solvable  for 
Icil  < e ]R^.  Finally  let  the  implication  (3.11),  (3.12) 

be  satisfied  with  ^^"^Aj^y^)  replaced  by  p)- 

Then,  for  h small  enough,  the  ^ from  (4.10)  satisfy 

(4.19)  = Aj^{y  + I h e^  ^^^(y)  + 0(h  ’^)  for  p < q 
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Hi 


resp. 


(U.20)  y + 0(h  ) /or  > q^. 


Proof:  For  h small  enough  , and  by  induction, 6 Bp(y). 

Analogously  to  (3.7)  and  (3.15)  we  have 


V Pn"l  V -V 

F(yj^_l>  = 1>  ^ F==(y^)  { I h 


(4.21) 


V q„  V -V  V” 

Po  -1  ^ ^ Po-1  ^ 

+ h*'^j;h  ^^e  Ay)  + 0(h  ^)  . 


Since  F”(y^)  exists,  the  problem 


P'(y„)v  = rcy^.j) 

is  uniquely  solvable  and  we  find, by  (4.21)  and  arguments  similar 


to  3 . , 


V Pn“l  “V 

v=h  i Z.  B e 

^=Po-i 


1 £-i^y)  + I ^ 

^ ^ i=p„ 


q„  V -V 


(4.22) 


+ 0(h  *') 


V 


,_i-y  + I B 


+ 0(h  "')  . 


Now  ('5’.^”(h,  ),F“(y  ))  is  W -admissible,  M is  linear  and 

h , o o 

e={E^,Ej°,  stable,  so  we  find  with  x := 

from  (4.10),  (4.17),  (4.18),  (4.21),  (4.22)  that 


^'h,P-r%,o 


"Pp-1  ''r'  "^'"Pp-1 

= T=h  Aj^{  I h ^x,£-l^y^  ^ 


'=pp-l 


+ I 

\=p„ 


e+  J + 0(h  '") 

\ , X. 
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o> 


>uk!  :.o  lin.Tlly 


'h,2. 


= 


= + I 


^=P( 


h ^ e ,(y)}  + 0(h  "') 

1 5 X, 


foi’  p , < and  (M.20)  for  . □ 

Instead  of  tlie  original  iteration  method  (M.7)  we  have  studied  in 
Theorem  4.3  tlie  method  defined  by  (4. 10). What  did  we  lose  by  this 
change? 

The  or  OP.  4.4:  Under  the  conditions  of  Theorem  4.3  one  finds^  using 
(4.7)  instead  of  (4.10),  again  the  relations  (4.19)  and  (4.20). 

In  (4.19)  only  the  e^  but  not  the  p^,qj^  are  to  be  changed. 

That  r.eans,  from  the  viewpoint  of  asymptotic  expansions:  The 
relatively  simple  method  (4.10)  works  as  well  as  the  relatively 
complicated  method  (4.7). 

Proot  : Since  we,  usually,  cannot  solve  £-1^  " ^h°^^^^£-l^ 

(4.7)  diirectly,  we  do  it  via  (4.8).  Let  us  substitute  ^) 

foi''  ^)  in  (4.8).  In  (4.8)  we  vjill  perform  only  a few  iterations 

So  tl'.e  following  argument  is  valid:  Comparing  (4.8)  and  (4.10)  we 

find  tdiat  = h,  ^ implies  o i fi'’om  (4.10).  Now  by  (4.23) 

o h,o^  ^1  h,£-l 


P£-^  V 


^£  V 


"h,t-r’^h,o  = ^h^  \ ^ I ^ 


^=Po- 


£-1 


^=Pc 


+ 0(h  ^) 
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To  compute  from  (M.lO)we  have  to  determine  ^rid  we  find 

with  . = t,  , 'I’v,(h>T  and  the  (v  ,v  )-smoothnc33  of 

Jl  i JlJljO 


*h«i> 


This  and  the  conditions  in  Definition  4.1  imply 


- ‘‘h";!’  - = ‘h  h'  h"’E  (y)}  * 0(h  "h 

'=P|l 


So  again  using  (4.G)  and  the  dSt  -admissibility  of  ^ ^ ) j 1 " (y^ ) ) 

we  find 


^ ^ ^ Qo 

(4.25)  1,^  - - - I h^g  .(y)+0(h 

where  the  g^^j^(y)  ^ \=p^(l)q^,  unless  tlie  i =p^,_^  (1  )pjj^-l  , 

are  reproduced  more  exactly  than  we  i^equired  and  unless 
additional  conditions  on  the  higher  derivatives  ai'^e  imposed.  If  we 
now  define  an  improved  n'ji^  analogously  to  (4.7)  and  (4.10)  by 

’^h,£-l  ■ ^^2"^h,o^  ^ ’^h,«-l  " 2"^1  ^ I’^^h  ,o^  ^’^h  , 2"^  1 ^ 

we  find  with  (4.24  ),  (4  .25  ) and  f,  p(v)  e . (y) 


V 


,,,  = ty  - I h ^e;=^^(y)} 


+ 0(h  ^)  . 


^=P( 
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and  further 


So  rij^"j^has  the  same  asymptotic  expansion  like  ^ 
iterations  via  (4.8)  would  give  the  same  result.  n 

Again,  like  in  3 . , it  is  not  necessary  to  know  F(y^_^)  in  (4.10) 
exactly, but  only  approximatively : If  we  use  in  (4.10)  an  approximation 
i(;''  to  F(yjj^_^)  analogous  to  (3. 17), we  find 


(nJ'  ) = 0,  so  n,"  = ^ > 

h h ,o  h ,o  h ,o 


(4.26)^  ^h(T^h;o)(^hU-l  ■ ^h'o^  = 


" ^h  + I ^ 'g,,^.l(y)>+0(h  J') 

^=P£ 


Ph,£  •"  Ph,£-1  " ^^h\'£-l  " '^h,o^  ‘ 


Theorem  4.5:  Define  hy  (4.26)  and  let  the  conditions  of 

Theorem  4.3  be  satisfied  with  £-i»  (3.11)  resp.  (3.12) 

replaced  hy  Then  (4.19)  resp.  (4.20)  are  valid  with 

Ph,£’  S,£  ^^Vlaced  by 

We  now  proceed  to  study  the  relations  between  the  results  in  3. and 
4..  (3. 9)  and  (4.10)look  very  alike  and  so  the  following  Theorem  is 

I 

j not  astonishing: 
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Theorem  4.6:  Let  be  a linear  discretization  method  and  let 

V . 

F'Ky^)  in  (3.9)  satisfy  (3.5)  and  ^ ^ 

with  the  ‘^y,"(hj^  q)  in  (4.10)  resp.  (4.26)  and  let  i|)  (Aj^u)  =i|i"  ( Aj^u ) + 0 (h  " "^ ) , 

If  we  write  the  results  of  (3.9)  resp.  (3.16)  as  riv.  a resp.  nr  « and 

n ) )6  n 5 it  j 

the  results  of  (4.10)  resp.  (4.26)  as  n^  ^ resp.  nj,'  then  j 


(4.27) 


^ ’ £ = 0,1,...  . 


Proof : With  the  denominations  in  Theorem  4.6  we  have  by  (4.10) 
resp.  (4.26) 


’^h,£  ” '^h,£-l 


^^^h,£-l  “ ’^h,o^  ’ 


^h,£-l  ^h  ,£-l ,o  ” ^h,o 


’^h,£  “ %,£-! 


^^h,£-l  ■ '^h,o^ 


and  therefore  by  (4.10)  and  (}),  (F"(y  ))  = ^'l:(n^,  ) + 0(h 

’h  o h 'h,o 


So, if  (4.27)  is  already  proved  for  £-l,w’e  have 


4'h^r"(yo^^^'^h,£  " '^h,£-l^  ""  '^li^^y£-l^  ^ ^ ^ > 
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and  so 


"h.H  = "h.l  * ’ • 


The  results  for  ri^  n obtained  in  the  same  way.  □ 

ri ) X/ 


The  idea  of  using  N.P.  goes  back  to  ZADUNAISKY  [13, 14, IS],  again 
discussed  by  STETTER  Cl2].  They  treated  initial  value  problems  of 
ordinary  differential  equations  where  STETTER  used  our  first  method, 
too.  The  method  of  N.P.  was  applied  by  FRANK,  HERTLING,  UEBERHUBER 
[4,5,6]  to  initial  and  boundary  value  problems  of  ordinary  diffe- 
rential equations.  In  tlie  next  papers  we  will  give  corresponding 
results  for  non-smooth  starting  values  y^  and  present  some  examples 
of  our  general  theor'y. 
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